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ABSTRACT 

The asteroseismic analysis of white dwarfs allows us to peer below their photospheres and 
determine their internal structure. At ~ 28, 000 K EC20058-5234 is the hottest known pulsat- 
ing helium atmosphere white dwarf. As such, it constitutes an important link in the evolution 
of white dwarfs down the cooling track. It is also astrophysically interesting because it is at 
a temperature where white dwarfs are expected to cool mainly through the emission of plas- 
mon neutrinos. In the present work, we perform an asteroseismic analysis of EC20058-5234 
and place the results in the context of stellar evolution and time dependent diffusion calcula- 
tions. We use a parallel genetic algorithm complemented with targeted grid searches to find 
the models that fit the observed periods best. Comparing our results with similar modeling of 
EC20058-5234's cooler cousin CBS 1 14, we find a helium envelope thickness consistent with 
time dependent diffusion calculations and obtain a precise mode identification for EC20058- 
5234. 
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1 ASTROPHYSICAL CONTEXT 

, White dwarfs are the end product of the evolution of around 98% 
' of the stars. Buried in their interiors are the records of physi- 
cal processes that take place during earlier stages in the life of 
, the star. Nuclear reaction rates during the core helium burning 

■ phase set the core composition of white dwarfs, while the relative 
time spent burning hydrogen and helium during the AGB phase 

. and massloss episodes determine the thickness of the helium layer 

■ jLawlor & MacDonaldl2006l : lAlthaus et al J2005h . 

Helium atmosphere white dwarfs (DBs) comprise roughly 
20% of the population of field white dwarfs, with most of the 
remaining 80% consisting of their hydrogen atmosphere (DA) 
cousins. The majority of white dwarfs in both of these spectral 
classes are thought to arise from the evolution of isolated main- 
sequence stars with masses that are insufficient to ignite carbon 
fusion, which ultimately leave their hot carbon/oxygen cores to de- 
scend the white dwarf cooling track. The bifurcation into two spec- 
tral classes is thought to occur during post-asymptotic-giant-branch 
(post- AGB) evolution when, in some cases, a very late thermal 
pulse bums off the residual hydrog en in the envel ope, producing 
a nearly pure helium atmosphere (Iben et al]|l983h . Such objects 
are then supposed to return to the white dwarf cooling track as PG 
1159 stars, which are widely believed to be the precursors of most 
DB white dwarfs. 

If we assume that there is an evolutionary connection be- 
tween the PG 1159 stars and the cooler DB white dwarfs, we 
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can look to several independent groups who have used time- 
dependent diusion calculations to follow the changes in the interior 
structure of these objects as they cool i Dehner & KawaleJ 1 19951 : 
iFontaine & Brassardll2002l : iMiller Bertolami et al.ll2006h . The hot 
PG 1159 stars, having recently emerged from the born-again 
phase, contain envelopes with a nearly uniform mixture of he- 
lium (He), carbon (C), and oxygen (O) o ut to the photosphere 
( iDreizler &HebeJl 19981 : iHerwig et al.ll 19991) . As they cool, the he- 
lium diffuses upward and gradually accumulates to form a chem- 
ically pure surface layer. Through the DBV instability strip, this 
process is still ongoing so that instead of a pure helium layer sur- 
rounding the carbon and oxygen core, one has a region where the 
carbon and helium are still mixed. This leads to a double- layered 
structure, with the pure He surface layer overlying the remainder 
of the uniform He/C/0 envelope, all above the degenerate C/0 core 
(see Fig[T). 

A key prediction of the diffusion models is that, for a given 
stellar mass, the pure He surface layer will steadily grow thicker 
as the DB star cools. The only available observational tests of 
this prediction come from asteroseismology - the study of the in- 
ternal structure of stars through the interpretation of their pulsa- 
tion periods. Helium atmosphere white dwarfs (DBs) are found to 
pulsate at effective temperatures r anging between 21,000 K and 
28,000 K iBeauchamp et al.lll999l : rCastanheira e t al.' 20051. There 
are currently 20 kno wn pulsating DBs (DBVs) (Nittaet al. 200^; 
iKilkennv et aLll2009l) . To date two of them, GD358 and CBS114, 
have been the object of detailed asteroseismic analyses. When it 
is well behaved, GD358 is the poster child for white dwarf astero- 
seismology, offering for analysis 1 1 dipole modes with consecutive 
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radial overtones k ranging between 8 and 18 JSradlev & Winged 
Il994h . However, an ambiguity between structure in the core and 
structure in the envelope has made it difficult to uniquely determine 
its core structure CMontgomerv et al. 2003). The observed pulsation 
spectrum o f CBS114 also c ontains 1 1 dipole modes, ranging from 
k=8 to 20. iMetcalfd ( |2005|) performed a systematic asteroseismic 
analysis using a genetic algorithm to search for best fit models of 
this star The models used in that work included realistic carbon 
and oxygen core abundance profiles that were allowed to vary to 
find the best fit. 

While the ambiguity in the interior structure of GD358 made 
it impo ssible to get a cl ear picture of its core chemical abundance 
profile. iMetcalf e ( 2007) was able to compare asteroseismic analy- 
ses of GD358, CBS 1 14 and EC20058-5234 that assumed pure car- 
bon cores. The analysis showed a correlation between the thickness 
of the helium layer and the effective temperature of the models, 
with EC20058-5234 having the lowest helium layer mass. While 
the result was consistent with the predictions of time dependent 
diffusion calculations, there is no physical justification to assume 
pure carbon cores. It is not consistent with stellar evolution calcula- 
tions and core structure influences asteroseismic fits. In the present 
study we perform a more physical asteroseismi c fit of EC20058 - 
5234 to compare with the CBS 114 fits done bv lMetcalfel ( |2005|) . 
where non-zero core oxygen abundances were considered. 

EC20058-5234 was the target of a Whole Earth Telescope ob- 
serving campaign (WET; Nather 1989) that revealed 8 or 9 stable, 
independent modes (Sullivan et al. 2008) . With that many modes, 
we can expect the asteroseismology of EC20058-5234 to yield use- 
ful constraints on its properties. We perform the first detailed aster- 
oseismic analysis of EC20058-5234 to determine its stellar param- 
eters and internal properties. In Sect. [2l we summarize the clues 
we used from spectroscopy to guide our asteroseismic analysis. In 
Sect. [3] we describe our models and the method we followed. We 
present our results in Sect. |4] and in Sect. |5] discuss our results in 
the framework of stellar evolution. We conclude in Sect.|6] 



2 CLUES FROM OBSERVATIONS 

2.1 Spectroscopy 

Assuming a pure helium atmosphere, iBeauchamp et all h999l) 
determined an effective temperature for EC20058 of 28,400 ± 
1,500 K and a logg of 7.86 ± 0.10. DB white dwarfs present a 
unique challenge to spectroscopists because of the lack of hydro- 
gen lines in their spectra. Because no spectrum is free of noise, it 
is possible to add in the models trace amounts of hydrogen that do 
not result in detectable lines in the spectrum. Unfortunately, even 
trace amounts of hydrogen can lead to significant differences in the 
inferred effective temperature and logg. Introduci ng a maximum 
trace amount of hydrogen N(He)/N(II) = -3.5, IBeauchamp et al.l 
determined an effective temperature of 27,100 ± 1,500 K and a 
logg of 7.80 ± O.IO. The spectroscopic "box" for EC20058-5234 
is therefore (25,600 K < Tes < 29,900 K,7.70 < logg < 7.96). 
With our models, the log g constraint translates to a mass constraint 
0.46 Mo< M.< 0.60 Mq. EC20058-5234 is a hot DBV and has a 
lower than average mass for a white dwarf. 

2.2 Pulsation spectrum 

EC20058-5234 was first discovered to pulsate by iKoen et al] 
l ll995b . In a 1997 Whole Earth Telescope run on this star. 



Table 1. Observed periods in EC20058-5234 



Mode name Period [s] Notes 



fl 


539.8 




f2 


525.4 




f3 


350.6 




f4 


333.5 




f5 


(286.6) 


m=+l rotational split of f6 


f6 


281.0 


High amplitude, stable mode. 






photometrically identified as f = 1 


f7 


274.7 


Possible m=- 1 rotational split of f6 


f8 


256.9 


High amplitude, stable mode 


fl6 


(207.6) 


m=+l rotational split of f9 


f9 


204.6 




fll 


195.0 





ISullivan et al.l ( |2008|) found 1 1 fundamental modes, listed in Ta- 
ble[I] The two highest amplitude modes f6 (281.0s) and f8 (256.9s) 
are remarkably stable and attempts have b een made to use them to 
measure a cooling rate for EC20058-5234 jPalessio et ai]|2010l) . 

I Sullivan et al.l found two modes each split by 70 //Hz (f6 and 
f9). If we assume these modes are C=l modes (a reasonable as- 
sumption from asymptotic period spacing arguments and the rela- 
tively low log g found from spectroscopy), then the 70 fiHz splitting 
is consistent with a rotation period of the star of 2 hours. The am- 
plitude ratios of harmonics and combination peaks to parent modes 
supports the hypothesis that f6 is an ^ = 1 mode (Yeates 20od). 
Since f9 has an identical frequency split likely due t o rotation, it is 
reasonable to assume that it is an € = 1 mode as well. lSuUivan et all 
also suggest that perhaps f7 is the third member of the (f5,f6,f7) 
rotationally split {=1 triplet, invoking the existence of a 3 kG mag- 
netic field to account for the uneven frequency splitting. f7 could 
also be a mode of its own that happens to lie close to where the 
m=+I member of the (f5,f6) multiplet would be if it were present. 
Preliminary work ( Bischoff-Kim 2008) did not find that fits exclud- 
ing the 17 mode (e.g. invoking a magnetic field) were significantly 
better than fits that included the f7 mode as an independent mode. 
In the 6 parameter asteroseismic fits presented in this work, we 
consider both alternatives. 



3 THE MODELS 

For this work, we used a parallel genetic algorithm applied to white 
dwarf asteroseismology (M etcalfe & Charbonneau.2003) . comple- 
mented with targeted grid searches. The genetic algorithm is an 
efficient way to search vast areas of parameter space to find local 
minima. Grid searches around these areas allow us to refine our 
results and produce a picture of parameter space around these lo- 
cal minima. To compute all our models, we used the White Dwarf 
Evolution Code (WDEC). 

The WDEC evolves hot polytrope models from temperatures 
close to I00,000K down to the temperature of our choice. Mod- 
els in the temperature range of interest for the present study are 
thermally relaxed solutions to the stellar structure equations. Each 
model we compute for our grids is the result of such an evolution- 
ary sequence. 
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Table 2. Parameters varied in tlie fits 
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Figure 1. Oxygen and helium abundance profiles for the Salaris like (solid 
lines) and ramp (dashed lines) best fit models. 



3.1 Input physics 

The WDEC is described in detail in lLamb & van Homl ( ll975h and 
IWooj ( ll99(]|) . We used smoother core composition profiles and ex- 
perimented with the m ore complex p rofiles that result from stellar 
evolution calculations jSalaris et al.lll997) . We up dated the enve- 
lope e quation of state ta bles from those calcul ated bv lFontaine et al.l 
( 1 19771) to those given bv'Saum on et al. I (fl995h . We use OPAL opac- 
ities jlgle sias & Rogers 1996) and plasmon neutrino rates pub- 
lished bv lltoh et al.l (1996). The new envelope equation of state ta- 
bles and pl asmon neutrino r ates are input physics we updated since 
the work of lMetcalfd l l2005h . 

DBV's are younger than their cooler cousins the DAV's. Time 
dependent diffusion calculations show that at 24,000 K, a typi- 
cal temperature for a DBV, the carbon has not fully settled into 
the core of the star yet, and we expect double layered helium 
layers, as shown in Fig. [T] (the chemical profiles corre sponding 
to our best fit model for EC20058-5234). Following iMetcalfel 
( l2005l) , we adopted and parameterized this structure in our mod- 
els. log(Menv) marks the location of the base of the helium layer 
and log(MHi;) marks the location where the helium abundance rises 
to 1. log(Mi.„v) and log(MHe) are mass coordinates, defined as e.g. 
log(Me„v) = - log(l-M(r)/M,), where M(r) is the mass enclosed in 
radius r and M, is the stellar mass. We did not treat the helium abun- 
dance in the carbon/helium region as a free parameter, but adopted 
the values predicted from time dependent diffusion calculations in 
iDehner & Kawaleil [m^) . 

There are two parameters associated with the shapes of the 



Parameter 


Range 


Description 


7'eff 


25000 - 30000 K 


Effective temperature 


M. 


0.450 - 0.575 Mo 


Stellar mass 


log(Menv) 


-2.00 to -4.00 


Location of the base of 






the helium layer 


log(MHe) 


-5.00 to -7.00 


Location of the carbon/helium 






to pure helium transition zone 


qftn 


0.10-0.85 M. 


Location of the edge of 






the homogeneous C/Ocore 


x„ 


0.0- 1.0 


Central Oxygen abundance 



oxygen (and carbon) core composition profiles: the central oxygen 
abundance (X^), and the edge of the homogeneous carbon and oxy- 
gen core (qfn,). We show an example of a basic oxygen abundance 
profile in Fig.[T]along with a Salaris-like profile. We tried both the 
simple profiles and Salaris-like profiles for the carbon and oxygen 
abundances, varying the parameters X(,and qfn,. 



3.2 Global exploration of parameter space with the Genetic 
Algorithm 

To find a best fit model to the observed periods of EC20058-5234, 
we varied 6 parameters in the ranges listed in Table 2. The effec- 
tive temperature and mass range were decided based on the spec- 
troscopy. 

The goal in white dwarf asteroseismology is to find the white 
dwarf model whose periods best agree with the observed periods. 
The genetic algorithm is able to explore the more promising regions 
of parameter space, and find the neighborhood of a global mini- 
mum. This supposes the existence of a global minimum. With 8 or 
9 periods for 4 or 6 parameters, the problem is usually well con- 
strained and the genetic algorithm successful. For EC20058-5234, 
the genetic algorithm was able to find a unique, global minimum. 

To calibra te the models used here against those used by 
iMetcalfdJioO^ , we started with pure carbon core models (Xo = 0). 
For such fits, 4 parameters describe the models fully: the effec- 
tive temperature, the mass, and the 2 envelope structure parame- 
ters log(Menv) and log(MHc)- We know from stellar evolution that 
nuclear burning in the phases leading up to the white dwarf stage 
resu lts in cores made up not only of carbon, but also of oxygen 
(e.g. lAlthaus et al ]|201 0). The C/O abundance ratio depends on the 
C(a,y)0 reaction cross-section. We also know (Montgomerv et al] 
that asteroseismic fits are most sensitive to chemical transi- 
tions in the core. In light of these facts, more recent aster oseismic 
fits have considered varied core compositions (e.g. Metca lfe et al.l 
1200 ll : iBischofT-KimlbOOgl) . We performed such fits, varying the 2 
core parameters Xo and qfm in addition to the 4 parameters of the 
pure carbon models. We tried both basic core abundance profiles 
with an artificial, linear decrease in oxygen with radius ("ramp" 
profiles) and m ore ph ysically realistic profiles, based on the work 
bv lSalaris et'al] imt l. 

As mentioned in Sect. |2.2| we also have two different hypothe- 
ses concerning the number of independent modes we have at our 
disposal to fit. If one believes that the 274.7 s mode is the result 
of an uneven rotational split of the 281 s mode, then we have only 

8 independent modes. If one rejects that hypothesis, then there are 

9 independent modes to fit. We tried both possibilities. For clarity 
and future reference, we summarize the different kinds of models 
considered in table 3. 
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Table 3. A summary of the different models considered 



Name 


Description 


Parameters 




PureC 


Pure carbon core 


Teff, M., log(Menv 


), log(MHe) 


Ramp 


C and core 


Teff, M., log(Menv 


), log(MHe), Xo, qfa 


Salaris 


Salaris et al. (1997) 








C and core 


Teff, M., log(Menv 


), log(MHe), Xo, qfm 



Table 4. Best fit parameters and goodness of fit. BIC (Bayes Information 
Criterion) is defined in the text. 



8 periods 


Model 


reff[K],M.[Mo] 


log(Me„v 


), log(MHe), X„ 


qfm 


0"RMS 


BIC 


Pure C 

Ramp 

Salaris 


29600 0.530 
28950 0.520 
29200 0.525 


-3.46 
-2.16 
-3.45 


-6.26 

-6.02 0.99 
-6.10 0.78 


0.32 
0.32 


2.25 s 
1.87 s 
1.63 s 


9.25 
9.77 
8.81 




9 periods 


Model 


reff[K],M.[Mo] 


log(Me,„ 


), log(MH.), Xo 


qfm 


0"RMS 


BIC 


PureC 

Ramp 

Salaris 


29650 0.530 
29000 0.515 
29300 0.510 


-3.44 
-3.50 
-3.50 


-6.34 

-6.90 1.00 
-7.10 0.86 


0.31 
0.32 


2.60 s 
2.04 s 
2.06 s 


11.3 
11.3 
11.4 



4 RESULTS 

We present the best fit parameters for each class of models consid- 
ered along with a measure of the goodness of each fit in Table |4] 
We constrained the 204.6 s and the 281 s modes to be f = 1. We list 
the periods calculated for the best fit models in the appendix, along 
with their mode identification, cr rms is defined as 

I T,""'" (Peak - PobsY ^, ^ 

0"RMS = A , (1) 

where Uobs is the number of periods present in the pulsation spec- 
trum. The index labeled "BIC" (Bayes Information Criterion) is a 
parameter that measures an absolute quality of the fit, by taking 
into account dif fering numbers of data points and free parameters 
(e.g |Liddl3l2007l) . It penalizes fits involving a greater number of 
parameters relative to the number of data points. It is given by 

BIC = Uobs In(cr^s) + np^u- In(nobs), (2) 

where ripa,- is the number of free parameters in the fit. 

From Table|4]we see that in general, the 8 period fits are better 
than the 9 period fits (even correcting for the fact that it is eas- 
ier to fit 8 periods than 9) and among the different core profiles, 
the best fit is the Salaris-like model (with 8 periods). Most models 
point to a high effective temperature as expected from spectroscopy, 
and a mass around 0.52 Mg. The best fits also have a rich to pure 
oxygen core, and log(Menv) — 6 to -7. In Fig.|2]we show how we 
approach the 8-period Salaris-like best fit model by taking slices 
through the 6-parameter space at the location of the best fit model 
(i.e. for the effective temperature plot, we fixed M. to 0.525 Mq, 
log(Me„v) to -3.45, log(MHe) to -6.10, Xo to 0.78 and to 0.32 
while allowing the effective temperature to vary around the best fit 
value of 29200 K). We show the chemical composition profiles for 
that model in Fig.[T] 




Figure 2. ctrms as a function of each of the 6 parameters for the Salaris, 
8-period best fit model. The vertical axis in each graph is ctrms in seconds. 
The symbols represent the models calculated, while the lines are polynomial 
fits to guide the eye. 



5 DISCUSSION 

It is interesting and comforting to note that no matter what the 
model details are, the best fit parameters are consistent. Parame- 
ters that dictate the locations of transition zones, in particular in the 
core (qfm) are especially well determined. The modes are sensitive 
to the location of chemical transition zones and to the sharpness 
of the transitions. The sensitivity of the asteroseismic fits to qf,^ is 
visible in the last panel of Fig.|2] Even with steps of 0.02, the mini- 
mum is still barely resolved. The fact that the best fits prefer a high 
oxygen central abundance (up to pure oxygen) is a sign that the 
modes are expecting a sharp transition at qf„, = 0.31 to 0.34. With 
our current parameterization of the core composition profiles, the 
only way to achieve sharper transitions is to increase the oxygen 
abundance as that allows it to drop more sharply down to zero. It 
could also explain the preference for the Salaris-like core profiles, 
as the drop in oxygen at qf,^ is inherently sharper in these mod- 
els (see Fig.[T]l. From our results, we can therefore conclude that a 
sharp composition transition happens at qf^ = 0.31 to 0.34, but we 
have a much weaker constraint on the actual central abundance of 
oxygen. 

In stellar evolution calculations, the central oxygen abundance 
for a given mass model is sensitive to the cross section adopted for 
the '^C(q', y)"© reaction (as well as the prescription for convec- 
tion). The '^C(a, 7)"0 reaction rates measured in the l ab also turn 
out to have fairly large uncertainties. iMetcalfe et alj ( l2002h find 
that these uncertainties propagate into uncertainties in the central 
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Table 5. EC20058-5234's best fit parameters in the context of testing time- 
dependent diffusion 



Parameter 



Pure carbon EC20058 models 
Metcalfe (2007) This work 



Salaris like models 
EC20058 CBS 114 



reff[K] 


28100 


29600 


29 200 


24900 


M.[Mo] 


0.550 


0.530 


0.525 


0.640 


log(Menv) 


-3.56 


-3.46 


-3.45 


-2.48 


log(MHe) 


-6.42 


-6.26 


-6.10 


-5.94 


x„ 






0.78 


0.71 


qfm 






0.32 


0.38 



oxygen abundance of ±0.1. There are also uncertainties in the cen- 
tral oxygen abundance asso ciated with the treatment of convection 
chosen (e.g lMetcalf j|2003 h. All uncertainties considered, the only 
strong constraint stellar evolution calculations place on the central 
composition of white dwarfs is that it is roughtly 50/50, with pos- 
sibly an excess of oxygen over carbon. With the standard reaction 
rates and treatment for convection, we expect the central abundance 
for a white dwarf of mass 0.530 Mp to be between ~ 0.7 and 
~ 0.9 jAlthaus et aT]|2010l : ISararis et al.lll997l) . Our present analy- 
sis is approximately consistent with these results, though again, we 
are more sensitive to the location and shape of the chemical transi- 
tion zone. 

Diffusion theory predicts that helium diffuses outward, as the 
heavier elements settle toward the center of the star under the influ- 
ence of gravity. With time, this leads to a thicker and thicker pure 
helium layer that sits on top of a mixed core (i.e. log(MHc) gets 
larger with decreasing effective temperature as the star cools with 
time). The main question we wanted to answer in this paper, was 
"Does EC20058-5234 being hotter and younger, have a thinner 
pure Helium layer than CBS 114?". The answer is yes (table [Sj. 
In reaching that conclusion, we must keep in mind that the models 
used for EC20058-5234 include updated physics. Can we compare 
the two? We performed 4 parameter fits so we could compare our 
results with previous asteroseismic fits of EC20058-5234 that were 
obtained with the same mod els that were use d for the 6 parame- 
ter CBS114 fits performed bv^ Metcalfj l l2005b . We place these re- 
sults side by side in Table |5] as well. While the parameters are not 
identical, we recover the same trends (of particular importance, the 
helium layer parameters are consistent). 

EC20058-5234 also has a lower mass than CBS 114 and 
we need to c onsider the stellar mass dependence of log(Menv). 
lAlthausetallllolO ) find that lower mass models have a thicker he- 
lium envelope. Just based on mass, this would mean that EC20058- 
5234 should have a higher (less negative) log(Me„v) than CBS 114. 
That is not the case and we can attribute the thinner helium enve- 
lope in EC20058-5234 to the fact that it is younger. 



6 CONCLUSIONS 

This study reinforces the fact that white dwarf asteroseismology 
is sensitive to core chemical abundances. In parameterizing com- 
position profiles, questions arise as to how many parameters are 
needed, how many we can afford to vary (the more we have, the 
more computationally intensive modeling becomes), and what are 
the best ways to parameterize the profiles. Because of these difficul- 
ties, it is tempting to ignore the importance of properly modeling 
the core chemical profiles, but the modes carry information about 
the core structure and we must decipher that information. Stellar 



evolution calculations are a good starting point and o ne reason we 
tried c hemical profiles that mimic the ones found bv lSalaris et al.l 
( Il997t) . More recently. I Althaus et al.l (12010') have performed calcu- 
lations that evolve stars from the Zero Age Main Sequence to the 
white dwarf stage, carefully treating massloss and time dependent 
diffusion of the elements. Pulsating white dwarfs that have over 
half a dozen modes like EC20058-5234 are good candidates to test 
these chemical profiles. 

In this paper, we used a thorough asteroseismic analysis of 
EC20058-5234 and compared it with a similar study (varying core 
parameters) of CBS 114 to test diffusion theory in white dwarfs. 
EC20058-5234 is also astrophysically interesting because it is ex- 
pected to be hot enough to lose a significant amount of energy 
through the emission of plasmon neutrinos. While the neutrinos 
are not detectable by direct methods, they have a significant ef- 
fect on the rate of cooling of the star. In turn, the rate of cooling 
can be measured because of its effect on the pulsations. The pe- 
riods grow longer over time as a result of cooling. The effect is 
very small (~ lO"'"* seconds per second), but detectable over a long 
enough time. Ideally one can measure an evolutionary cooling rate, 
as has been done for another pulsating white dwarf, G117-B15A 
(Kepler etal. 2005). Recently, pulsational data collected over a 
time period of 13 years for EC20058-5234 was assembled and an- 
alyzed to attempt to measure the coolin g rate of that star, b ut the 
measurement turned out not to be trivial i Dalessio et al j2010!).One 
hypothesis about the strange results is that rotation may play a sig- 
nificant role in the pulsations of the star. This would be consistent 
with the fact that EC20058-5234 appears to be a fairly fast rota- 
tor, with a 2 hr period. Typical (non-magnetic) white dwarfs are 
known from aste roseismology to rotate with periods ~ 1 day (e.g. 
iKepler et al .''2003) . The precise mode identification and improved 
interior models presented here may help shed light on the strange 
time evolution of EC20058-5234's pulsation periods. 
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Table Al. Periods and mode identification for the 8-period fits 
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279.4 
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349.1 
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256.9 
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lie A2. Periods and mode identification for 9 period fits 
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APPENDIX A: MODEL PERIODS AND MODE 
IDENTIFICATION FOR EC20058-5234 

We list here the periods calculated for the best fit models presented 
in the main text, along with their mode identification. The informa- 
tion contained in this appendix shows in more detail the quality of 
the fits and may prove useful for further asteroseismic studies of 
EC20058-5234. We stress again that the 204.6 s and 281.0 s modes 
were constrained to be f = 1 while the others were allowed to be 
either ^ = I or ^ = 2. 

For all best fit models, the 350.6 second mode has a good fit 
either as f = 1 (k = 5) or f = 2 (k = 13). In all cases, the two periods 
are less than 2 seconds apart and in most cases, less than a second 
apart. In other words, the two identifications are not distinguishable 
from one another. We chose to list both possibilities in tables lAll 
and lA2l (so the 350.6 second mode appears twice in each table even 
though it is a single mode). 



